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Abstract: Next-to- leading order corrections are an important aspect in the extraction 
of the Cabibbo- Kobayashi-Maskawa matrix elements |V^b| and |Vub| at i?-factory exper- 
iments: virtual and real photons couple to all charged particles of the decay, alter the 
resulting decay dynamics and enhance the total decay rate. We present a phenomeno- 
logically motivated model and a Monte-Carlo implementation to study electromagnetic 
radiative corrections to differential and total decay rates for semileptonic S-meson decays 
into exclusive pseudoscalar and scalar final states and apply it to B ^ Dlv and B ^ nlv, 
and B — )• Dq I v decays, respectively. We study such corrections with a phenomenological 
picture of point-like mesons (with some structure-dependent corrections), which is valid in 
the low-energy region of the photons and whose results we extrapolate over the complete 
phase space. The largest quantifiable uncertainty is due to the approximative matching to 
the Standard Model. In addition unknown structure-dependent contributions and model 
dependencies might have an impact on our findings. 
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1. Introduction 

In the Standard Model, the Cabibbo-Kobayashi-Maskawa (CKM) matrix Q governs the 
weak transitions between the up- and down-type quark generations. The precision deter- 
mination of its matrix elements and CP violating phase is the focus of intense research 
over the past decade. The combination of such measurements to test unitarity of the CKM 
matrix as a whole is considered a strong instrument for the search of physics beyond the 
Standard Model |, |. 

In this paper, we present a calculation and a simulation tool for the electromagnetic 
one-loop corrections in exclusive semileptonic S-meson decays based on a phenomenological 
model, which is applied to B ^ Dlu^ B ^ n I u, and B — t- Dq I v decays. Next-to-leading 
order corrections to such decays are an important aspect in the extraction of the CKM 
matrix elements |14b| and | V^bl at Z?-factory experiments: virtual and real photons couple to 
all charged particles of the decay, alter the resulting decay dynamics and enhance the weak 
decay rate. To correct for the altered decay dynamics, experimentalists use approximative 
all-purpose next-to-leading order algorithms. These exploit the fact that such corrections 
become universal in the soft photon energy limit [Q, ^. In addition, the total leading order 
decay rate of semileptonic decays are corrected by the known leading logarithm of the 
virtual corrections at parton level |7|. As implied by Q such a treatment might not be 
sufficient: corrections beyond the very soft limit alter the leading-order decay dynamics and 
can potentially alter acceptance studies of experimental results - and therefore measured 
values for |V^b| and |Vub|- 

Experience from exclusive semileptonic X-meson decays illustrate the importance of 
having a good understanding of radiative effects: until 2004 the global average of the 
extracted value of |Vus| from K^^ and -fC^g decays implied the violation of CKM unitar- 
ity by two standard deviations Q. Further measurements proved dissonant with these 
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findings |1C, 11, 12, 13 1, indicating that the achieved experimental precision needed an 
understanding of electromagnetic corrections at the percent level. Since for many decays 
next-to-leading order calculations do not exist, experiments often use the approximative 
all-purpose algorithm PHOTOS to study the reconstruction efficiency and acceptance. 

The accuracy of this approach was tested by the KTeV collaboration, using the measured 
photon spectra from radiative K^^ decays: the angular distribution of the simulated pho- 
tons did not agree well with the predicted spectrum [|l^. This lead to the development 



of the next-to-leading order Monte Carlo generator KLOR (see [16|), whose next-to-leading 
order calculation is based on a phenomenological model. Its predicted angular photon dis- 
tribution agreed satisfactorily with the measured spectra. Although using next-to-leading 
order calculations are more preferable than approximative all-purpose algorithms, they 
are often complicated to adopt for an experiment: electromagnetic next-to- leading order 
calculations only exist for a few decay modes, sometimes only valid in a limited region of 
phase-space. 

Over the last 10 years, an increasing amount of data and a better understanding of 
detector effects leads to a very accurate picture of physics at both S-factory experiments 
BaBar and Belle. This increased precision lead to the demand of knowledge of next-to- 
leading order electromagnetic effects beyond the precision of approximative all-purpose 
algorithms. Our paper was written to help closing this important gap and improving the 
status quo in a twofold way: first by providing a Monte Carlo generator that produces 
differential decay rates based on a next-to-leading order calculation valid in the low-energy 
limit of phase space with respect to the photon; and second by predicting the unknown 
isospin breaking contributions to the next-to-leading order enhancement of the total decay 
rate. The former can be used to conduct acceptance studies to improve the understanding 
of the role of such radiative events in the extraction of CKM matrix elements. The latter 
result give small corrections to the known parton level result, and therefore to the extracted 
values of \Vch\ and \V-ah\ from semileptonic decays. The uncertainties due to corrections 
beyond the soft limit for the latter are hard to estimate, and we restrict ourselves to conduct 
studies within the phenomenological model itself. The error budget to the number should 
therefore be used with care. 



Our model is based on the considerations of |16] and makes use of some of the formal 
ideas of |17|: The latter calculates the real next-to- leading order effects of semileptonic Ki^, 
decays. We use the formalism therein to split the radiative hadronic current in contributions 
that depend on the tree-level hadronic current and further corrections. The virtual one-loop 
diagrams are determined in a phenomeonological model are matched to the parton level 
correction. Our paper proceeds as follows: Section ^ briefly reviews exclusive B ^ Xlv 
decays at tree-level. Section ^ thereafter develops the next-to-leading order formalism 
and introduces the phenomenological model. Section |^ present the next-to-leading order 
matrix elements and the next-to-leading order differential decay rate. Section ^ describes 
the numerical methods used to integrate the studied next-to-leading order matrix elements. 
Section |6| presents the results for the next-to-leading order total and differential decay rate 
predictions, and Section ^ summariz es our result and presents our conclusions^. 

^Charge-conjugated modes are implied tiiroughout this paper and that natural units, c = 1 and h — 1, 
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2. Phenomenological tree-level decay revisited 





Figure 1: The tree-level weak B Xlv decay using Fermi theory is shown: the left hand side 
shows the phenomenological picture of mesons; the right hand side the decay at parton level ignoring 
the light spectator quark. The shaded circle in the transition represents the effective Fermi vertex. 

In the Standard Model the tree-level weak b — )• xlu decay in Feynman gauge is described 
by the transition matrix element^ 



= -i^V,^[u{p.)^^PLvipi)]h^ + o(-^) , (2.1) 

with the left-hand projection operator Pi^ = 1— 75 and the quark current /i^ = [x{q') 7^ Pl b{q)] . 
Since the momentum transfer \/t through the W propagator is very small compared to the 
T^-boson mass mw, terms of O (t/m^r) can safely be neglected and the full amplitude can 
be approximated by the first term in Eqn. (2.1). This term can also be obtained by an 
effective Hamiltonian 

71 



Ties. = Kb [i^u 7^ Pl i'l] K ' (2-2) 



with the quark current hn in position space. At leading order in the W propagator expan- 



sion the effective Hamiltonian Eqn. (2^) and the Standard Model yield identical expres- 
sions for the b ^ xlv transition amplitude. In order to describe the weak decay of the 
bound state of a P-meson, the quark current in Eqn. (|2.2| ) needs to be replaced with a 
hadronic current describing the hadronic matrix element. The resulting matrix element 
can be described by means of an effective Hamiltonian containing the phenomenological 
meson fields for the initial i?-meson and the final state X-meson. For the B ^ Xlv decay 
into a (pseudo) scalar final state the effective Hamiltonian can be written as 

T^eff. = ^ Kb i^u 7m Pl i^i [(/+ + f-)4>d^^B + (/+ - f^)<pB 5^</>] + • • • , (2.3) 

where the ellipses denote further higher order operators. The relevant hadronic degrees of 
freedom are the initial state P-meson and the final state X-meson and can be described 



are used. 

^In the following, A^m i^m) denotes the matrix element (total decay rate) at 0{a^) with m photons in 
the final state, the tree-level matrix element (decay rate) consequently is indicated as Mq (fJ!) . 
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by two scalar fields (pB and (j), respectively. The resulting amplitude for a i? — t- Xlu decay 
from the leading term in Eqn. ( p.3| ) is 

Ml = -i^ Hpu) 7^^ Pi. v{pi)] (2.4) 

with 

= {X{p)\ \B{pb)) = {pB+p)^f+ + {PB-P)^f- . (2.5) 



The effective Hamiltonian Eqn. ( |2.3| ) reproduces the correct tree-level matrix element. In 
the hadronic current Eqn. (2.5) the form factors f± = f±{t) are functions of the four- 



momentum transfer squared of the i?-meson to the X-meson system, e.g. 

t=[pB-p)\ (2.6) 

and describe the influence of the weak decay on the bound state quarks caused by the strong 
interaction. This dependence can be made explicit in Eqn. ( p.3p by adding appropriate 
higher order operators. As a consequence the effective Hamiltonian Eqn. ( ^.31) cannot 
teach us anything about the consitutents of these mesons, but it proofs useful to study 
low-energy electromagnetic corrections that do not resolve the bound state structure of the 
involved mesons. The tree-level differential decay rate in the i?-meson rest frame is given 

by 

dr{] = —4 \M^o\^dEdEu (2.7) 

with E = p^ and Ei = p^, and the explicit expressions of f± as a function of the momentum 
transfer squared ioi B ^ Dlv, B ^ Dq I v, and B ^ irlv decays can be found in App. 



3. Next-to-leading order formalism and phenomenological model 

The aim of this section is to develop a formalism to calculate the corrections at O^aGp) 
for B ^ Xlv decays and introduce a phenomenological model to perform the actual calcu- 
lation in an approximative manner. First the formalism used for virtual corrections and its 
connection to renormalization is discussed. Thereafter the splitting of the electromagnetic 
coupling into parts that can be described by minimal coupling to the leading term in the 
effective Hamiltonian Eqn. ( p.3| ), and terms that arise from further (known and unknown) 
higher order operators, is discussed. In the end of this section our strategy to obtain the 
approximative corrections at 0{aG-p) is summarized. 

3.1 Formalism for virtual corrections: short- and long-distance contributions 

Electromagnetic loop corrections to the B ^ X Iv decay amplitude modify the Hamilto- 
nian Eqn. ( |2.3D by adding further terms of 0{aG-p)- Such corrections to the Hamiltonian 
are calculated by comparing amplitudes in the full theory with the VK-boson present as a 
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dynamical degree of freedom to amplitudes in the effective theory with the W^-boson re- 
moved as such. The arising corrections come from regions of loop momenta of order ^ mw, 
since the effective Hamiltonian has been constructed in a way to correctly reproduce the full 
theory for momenta much smaller than mw- The standard approach to proceed involves 
calculating the one-loop graphs for the B ^ Xlv decay in the effective theory with counter 
terms and compare it to the renormalized Standard Model result. Fixing the counter-terms 
results in the desired matching of both results. The effective theory itself is nonrenormal- 
izable, but the Standard Model can be renormalized to measured quantities, e.g. Fermi's 
constant, in order to produce finite predictions. The arising infrared divergencies in such 
one-loop graphs are canceled with divergencies that show up in the infrared limit of of the 
B —7- X luj decay. Such a matching procedure was carried out in much detail by |jl^ for 
semileptonic Kaon decays. For semileptonic i?-meson decays usually the result of Q is 
used to correct in a similar manner for virtual corrections of 0{a Gp) of the quark current 
h^. Corrections that do not resolve the involved hadrons to parton level are small, but in 
order to study the full real and virtual corrections at 0{a Gp) consistently they need to be 
included. 

In order to proceed we assume that further (unknown) higher-order operators in Eqn. 
(|2.3| ) can be constructed in such a way, such that the full Hamiltonian reproduces the 
virtual electromagnetic corrections at all energies - explicitly that its high energy limit 
reproduces the parton level correction to Eqn. (|2.2D , i.e. the findings of ||6|, 0]. The ellipses 
in Eqn. ( |2.3D then contain two classes of higher order operators: derivatives of fields that 
derive the tree-level t dependence, and operators that describe the corrections at 0{a G-p) 
and maintain the transition from the phenomenological degrees of freedom of hadrons to 
the partonic degrees of freedom of quarks at high-energies, i.e., the Standard Model picture. 

The large separation of scales of the masses of the involved decay products in contrast 
to the ly-boson mass, explicit mx^fnB ^ ttiw^ allows for an alternative prescription to 
match such an effective field theory to the Standard Model and to obtain a meaningful 
result with respect to renormalization. A general divergent A^-point tensor integral of 
order p in the effective theory has the form of 

, fe^i . . . k^'-p 



T^^^-^^{p^,...,PM-i)^ dk'—-^ , (3.1) 

J do ■ ■ ■ dN-i 

with denominators do = k'^ and di = [{pi — k)'^ — mf^. This integral can be split by an 
Euclidian regulator /io ^ mx ,mB <^ mw into two regions 

T'''-''^{pi,...,PN^i) ^ dk^- +/ dk^- , (3.2) 

Jx2 do ■ ■ ■ dN~i do ■ ■ ■ dN-i 

where we introduced an infrared regulator A. Choosing the Euclidian regulator near mx 
and mB ensures that in the region k'^ < the relevant degrees of freedom are the bound- 
state mesons, i.e., the arising dynamic can be described by minimal coupling to the leading 
term of Eqn. { \2.3\j and the tree-level operators. Large photon momenta resolve the bound 
state quarks and the corrections from operators that restore the partonic degrees of free- 
dom will no longer be negligible. This qualitative argument is true up to intermediate 
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resonances, that can show up in the energy regime of k'^ < ^g, and breaks down at loop 
momenta ^ ij,q, when the virtual photon starts to resolve the bound state meson. Note in 
particular that the first term of Eqn. ( |3.2D contains only infrared divergencies, whereas the 
second integral is ultraviolet divergent only. Using this splitting the virtual next-to-leading 
order matrix element decomposes as 

-^0 = -^o,id(/"o)+A^o,sd(/^o), (3.3) 

where A4q i^{fJ.o) and Mq sd(^o) denote the contributions from the first and second integral 
in Eqn. (|3.2| ), respectively. Introducing an Euclidian regulator will in general destroy gauge 
invariance of the electromagnetic interaction. A splitting that leaves gauge invariance intact 
is given by 

rpl^l-t^P^p p ) s; n dk^ ( "' "' \ 

Jx-2 ydo...dN-i {k"^ - ^l) di . . .dN-i j 
Jxi [k^ - )U5j di . . . dN~i 



The second term in the first integral of in Eqn. ( |3.4D acts as an additional convergence 
factor that leaves the low-energy behavior intact, but suppresses the exchange of virtual 
photons of momentum larger than what introduces an unphysical photon- like vector 



field with opposite norm, as proposed by Pauli and Villars |19|. The photon mass //q 
leaves the high-energy behavior of the second integral of Eqn. ( |3.4D intact, but cuts off the 
exchange of particles with momentum lower than ^o- This again renders the first integral 



of Eqn. (3.4) ultraviolet finite, and the second integral of Eqn. (|3.4D infrared finite. One 



interesting feature of the splitting of Eqn. ( p.4| ) lies in its properties under renormalization: 
The term A^J g(j(/^o) carries the full ultraviolet behavior of M.\ and thus can be used to 
renormalize all parameters. Consequently, because the splitting is exact, all parameters in 
TWj id(^o) ^-re renormalized automatically. In the following we will refer to contributions 



like the first integral of Eqn. (3^) as long-distance (Id) contribution. If the matching scale 
is chosen such that ^ niB or mx, the electromagnetic interaction at this energy cannot 
resolve the underlying structure of the involved mesons, and only feels the scalar long- 
distance coupling - up to corrections from excited intermediate states and the apparent 
break down at loop momenta close to /^o- Contributions like the second integral will be 
denoted as short- distance (sd) contribution to the amplitude, that can resolve the partonic 
structure of the involved mesons. 

3.2 Formalism for real corrections: inner bremsstrahlung and structure-dependent 
terms 

In this section we will clarify the idea of long- and short-distance interactions in a more 
formal way in order to make it applicable for real emissions by splitting the electromagnetic 
current of the hadronic system into two components: inner-hremsstrahlung (IB) contribu- 
tions, which account for photon radiation from the external charged particles (i.e. described 
by scalar QED) and are completely determined by the non-radiative process, e.g. the lead- 
ing operator of Eqn. ( p.3| ) and further tree-level operators; and structure- dependent (SD) 
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contributions, which describe intermediate hadronic states and represent new information 
with respect to the IB contributions, e.g. from additional higher order operators of Eqn. 
( |2.3| ). This approach will allow us to split the real corrections correspondingly into IB and 
SD contributions. 

Coupling the electromagnetic current to the semileptonic S-meson decay amplitude 
yields 

i e ^V^^u{p,) Pr 7M - {Y^ + 2pf) + V - V ) viPi) > (3-5) 



with the hadronic current H^, as introduced in Eqn. (2.5). The hadronic vector and axial 



form factors of the B / X ^ coupling are given by the unknown non-local operator 

V - V = / d'^e'"'" {Xip')\T [/i^(O) Jr(^)] Wb)) ■ (3.6) 

It is J^™ the electromagnetic current and the quark-current in position-space. Eqn. 
(|3.6|) obeys the electromagnetic Ward-identity [pO[: 



r^^, = 0. (3.7) 
Gauge invariance of the resulting amplitudes reveals a connection between the tree-level 



and the next-to-leading form factors: Low's theorem |21|, ^] states that the leading contri- 
butions to the next-to-leading order amplitude in powers of the photon four-momentum fc, 
explicit and , are completely determined by the on-shell form factors of the tree- level 
decay. The work of [17| presents a set of arguments to include corrections beyond 0{k^), 
by separating the non-local operator Eqn. (^) into SD and IB contributions, as depicted 
in Fig. ||. We briefly restate them: 

• In order to describe two different physical mechanisms, the IB and SD amplitudes 
must be separately gauge invariant; 

• The SD amplitude contains terms of order k and higher. 

The second condition does not prevent the IB amplitude from containing terms of order k 
and higher. Splitting the amplitude under this restriction allows us to collect more terms 
in the IB part, by still using only the knowledge of the non-radiative matrix element. This 
offers the advantage to obtain a more precise predictions for the decay process, without 
formulating the (mostly unknown) SD contributions. The splitting of the transition matrix 
element requires a corresponding splitting of the non-local operator Eqn. ( |3.6| ) into SD and 
IB parts. 

We first consider the axial contributions: They are strictly zero for photon emissions 
producing a (pseudo)scalar intermediate B- or X-meson state b) oi Fig. ^, and therefore 
can be considered purely SD. They can be written in the form [^3| 

A^^ = A^^ = -i Sf^upa {Ai pP k" + A2 kP{pi + p^Y) 

-iEuXpaP^kP {pi+PuY (^3 {pi+Pu)^ + Aip^ . (3.8) 
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Figure 2: The coupling of tlic electromagnetic current to the hadronic and leptonic currents is 
shown: a) structure-dependent , b) scalar coupling to a intermediate particle with a B- or X-meson 
pole, and c) the leptonic coupling. The shaded circle represents the one-particle irreducible graph 
of the weak decay. 



Note that the Lorentz-invariant scalars Ai are non-singular in the vanishing photon-energy 
limit by construction and are functions of the three independent scalar variables that can 
be constructed out of pB-, p, and k. 

The decomposition of the vector current is 

V = ^M^ + ^5?' (3-9) 
where the IB piece is chosen in such a way, that within Eqn. ( ^.71) 

k'^V^^ = H^, (3.10) 

and therefore for the SD contributions it is imperative that 

fe%l° = 0. (3.11) 



The decay amplitude separates now as 



+ i e ^V^^u{p,) Pn 1^ I V^^ - Af, ]v{pi) . (3.12) 



It remains to construct the explicit form of , what is done in Sec. 11 for the real 
emission process. For the B — t- IB contributions, one finds 

r y'^ = — ^ 7^ [pb], + [pb), 

^ Pb • k ^ pb • k ' ^ 



{2i) + mi) {pb)^ r {m) {pb)^ - /s Iv 



with 



PB ■ k ^ '"^ pB-k 

+ {2i> + mi) Zf^^ + {mi) Zf_^, (3.13) 



Zf^ u = ^^-^ - f±{t)) , Zf, = F Zf, . , 

• [PB -P) 

/3 = (/+(t) + /_(t)), (3.14) 
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Furthermore t' = (pB — p — A;)^ is the radiative four- momentum transfer squared. Similarly 
the X — )■ X'j IB contributions are 



r 



with 



p ■ k p ■ k' 

^ (2|iB - m) Pu ^ im) Pv + /2 iv 

p ■ K p ■ k 

+ {2;^B - m) Zf^^ + (mi) Zf_^ , (3.15) 



h = if+it) - f-{t)) . (3.16) 



The SD vector contributions can be written as |24] 

^/f ° = Vi (A:^ - p • /c 5^1.) + V2 (A;^ (p« + Pi.) ^ - A: • (p; + p,y) g^u) 

+ V'3 U- (pi+Pi.) {pi + Pu) ^Pu - P ■ k {pi + Pu) ^{pi + Pu) , 



+ V4{k- {pi+pu)pf,Pu-p-kpf,{pi+pu)J , (3.17) 

where the Lorentz-invariant scalars Vi are functions of the three independent scalar vari- 
ables that can be constructed out of pb,P, and k. The IB vector currents Eqns. ( 3.13| ) 



and ( |3.15 ) satisfy the Ward-identity Eqn. ( |3.1C ) and are constructed in such a way, that 



all infrared (IR) singular parts in the limit oi p ■ k = 0, or pB ■ k = 0, respectively, are 
contained within. 



An alternative choice for the IB and SD separation used in this work is given by [25 
based on |26]. It can be obtained by a change of basis in Eqn. ( p. 17 ), and correspondingly 
shifting terms of 0{k) and higher into the SD contributions^. 

3.3 Phenomenological model 

So fare we did not introduce any model dependence and assuming one would have full 
knowledge of the effective Hamiltonian, that governs the semileptonic transition at 0{a Gp), 
one could now go ahead and perform the necessary calculations following the baselines for 
the real and virtual corrections outlined in the last two sections. However, the knowl- 
edge of the operators that connect the short-distance and long-distance domains, i.e. the 
SD contributions, is modest: The recent work of addresses the real SD corrections to 
B — )• Dlv^ decays by using lattice results of the D* — )• D 7 coupling to model the first 
dominant SD contributions, in the region of phase-space, where the D* is an on-shell in- 
termediate resonance. The authors of [|2^ discuss the matter for i? — )• vr / 7 decays, using 
soft-collinear effective theory (SCET) to isolate the expressions for the SD contributions 
in the soft-pion and hard-photon part of phase-space. The SD corrections to i? — )• I v 7 
are unknown, but given the large decay widths of the D^- and L'J-mesons, a considerable 



'Note that there is a subtlety in the used notation: IB and SD will always refer to V^|^, and — 



what should not be confused with the particular choice of the IB terms, Eqns. (B.13) and (3.15), which 
incorporates a number of structure- dependent terms. 
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correction to the low-energy result of the pure IB prediction can be expected. This is 
a daunting starting point, but nonetheless allows one to examine the low-energy limit of 
the effective Hamiltonian, where SD contributions play only a minor role, and extrapolate 
these findings to the complete phase space. Note that this Ansatz - which we describe as 
'phenomenological model' in the following - is superior to the use of all-purpose algorithms, 
that make use of the universality of electromagnetic corrections in the soft-limit: Our con- 
siderations are based on a (at last in principle) full process dependent next-to-leading order 
calculation valid near the soft-limit and extrapolated to complete phase space. 



3.3.1 Virtual corrections 

In order to obtain finite predictions for the overall next-to-leading order rate, the second 



integral of Eqn. (3^) describing the short-distance interaction has to be regularized and 
renormalized to the Standard Model, all parameters describing the long-distance interaction 
are then renormalized automatically. This is exact as long as the complete Hamiltonian 
Eqn. with all higher order operators is used to calculate both the short-distance and 

the long-distance parts. In practice however, as discussed before, this is not possible since 
the higher order operators moderating the transition from long-distance to short-distance 
interactions are unknown. Hence supposing the matching scale /io is chosen in such a 
manner, that it effectively separates the long- and short-distance regimes, the long-distance 
corrections can be calculated using the leading operator of the effective Hamiltonian Eqn. 
( |2.3| ) and the known tree-level higher-order operators: in this phenomenological model 
the Standard Model and its effective degrees of freedom are described by scalar fields, 
what is justified at low scales. For the short-distance corrections A^J sd ^^^^ Standard 
Model has to be invoked. This mere fact, however, directly leads to inconsistencies at 
the matching scale /io, where both models should give the same answer. This renders 
the matching approximative and the impact of the choice of the matching scale can be 
studied by varying the matching scale /io- The long-distance virtual corrections are further 



discussed along with the renormalization of the short-distance correction in Section 4.3 



3.3.2 Real corrections 

In order to obtain predictions for the real next-to- leading order corrections, we make use 
of the previously introduced separation into IB and SD contributions. We will neglect the 
largely unknown SD contributions. This is justified in the low-energy region of phase-space, 
which is dominant for the emission process. Only at high energies the emitted real photons 
resolves the bound state nature of the involved mesons, and the SD contributions become 
dominant over the IB ones. This reasoning, however, excludes intermediate excited states. 
In order to study the impact of such, the pure IB result presented in this work is compared 



to the findings of |H] and |25| in Sec. 3.1.4. Note that neglecting any SD contributions 



results in a model dependent choice of the IB contributions, what is discussed along with 
the construction of Eqns. ( |3.13 ) and (3.15) in Sec. 4.1. 
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4. Next-to-leading order calculation 



In this section Eqns. ( 3.13| ) and ( |3.15| ) and all long-distance one-loop graphs are derived 



using the phenomenological model outlined in the previous section. In addition the approx- 
imative expression for the summed long- and short-distance next-to-leading order decay 
rate is obtained. 

4.1 Real next-to- leading order corrections 



PB •P^- 






Figure 3: The real emission contributions are shown: a) and b) emission from point-like mesons, 
c) from structure-dependent contributions and d) from the lepton leg. In all shown diagrams the 
shaded circle represents the Fermi coupling. 

The real emission from all external legs and structure-dependent contributions are 
depicted in Fig. ^. With the photon coupling to the lepton leg, as given by Eqn. (^), 
the corresponding emission matrix element d) is 

M!=.e r P. ( - + ^^f^) , (4.1) 

where e* = e* (k) denotes the polarization vector of the real final state photon. It is further 
discussed along with the infrared regulator A in Sec. 

The emissions from the hadronic leg now have to be constructed in such a way that 
Eqn. ( ^l0|) can be fulfilled. The intermediate propagator, a), with a four-momentum 
squared (ps — k)^ generates a pole at niB, which corresponds to an intermediate S-meson 
state. We isolate the contribution of this pole in the hadronic emission matrix element: 

mI =ie^V,MP^)j^PLv{pi) I P^-'* H^(p^-k,p,t')-%,e*A ,(4.2) 

V2 \ PB ■ k + ^X^ J 

with t' = (pb — p — k^"^ . Similarly the intermediate propagator, h), with a four-momentum 
squared {p + k)"^ generate a pole at m, corresponding to an intermediate X-meson state. 
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Isolating the pole in the hadronic emission matrix element results in 



mI =ie^V^^u{p,)rPj^v{pi) (^-^-^^^H^{ps,p + k,t')-V^,e*''j . (4.3) 

The form factors in the hadronic currents H^(pB,p + k, t') and Hfj_{pB — k,p, t') are a priori 
unrelated to the ones in the tree-level hadronic current H/^. They are 

H^,{PB - k,p, t') = [pB +p)J+{t') + [PB - p)^ f-{t') - hit') , (4.4) 
H^{PB,P + k, t') = [pB + p)J+{t') + [pB - p)J-{t') + k^ hit') . (4.5) 

with f± being the off-shell form factors. It is easy to see, that the off-shell form factors can 
only differ from the on-shell ones by contributions that are regular at k = 0, which can be 
reabsorbed in the SD terms within V^i,. This is true under the assumption that the limit 
limfe_^o /± is well defined and V^i, does not possess any poles in A; = . Then 

A;^ V'' = H'" + k^ V"'' = HI" , (4.6) 

with H'lj^ = H^ipB — k,p, t') or H^ipB.,p + k, t'). Since ki, V'^'^ is regular, it cannot contain 
any terms proportional to i?'* and it must be 

f±it') = f±it) + kaZl, (4.7) 

with some arbitrary functional which we are free to choose. We construct in such a 
way, that /±(t') = f±it') and the matrix element describes the emission diagrams a) or b). 
The freedom of choosing Z^^ is nothing more than deciding what contributions we consider 
as a self-contained emission diagram. Another choice would just move contributions from 
a) and b) into or vice versa. It remains to prove that V^,y can be constructed in such a 
way, that it possesses no poles. The Ward identity implies 

lim k^ V^"" = lim F'^ + hm A;^ V^"" = H^" . (4.8) 

If is regular at limk^Q ky V^^ = then our assumption is trivially true. We therefore 
assume, that V^" has a pole at /c = and write limjt^o k^ V^^ = V^: 

lim ky V^"" = ipB +p)'" lim Uit') + ipB -p)^ lim /_(*') + V^' , (4.9) 

fc^O fc^O fc^-O 

= (pB+p)^/+(t) + (pB-p)^/-(t), (4.10) 

where we made use of our premise, that the limit limjfc_^o / is well defined and finite. When 
limfe_^o /±(*') = f±{t) it follows = Q and our initial assumption is true. We therefore 
assume limfe_>o /±(t') / /±(i), then 

= iPB+pY [hit] - hit)) + ipB-pY {f-it) - f-it)) . (4.11) 
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We can split V^'^ into its singular and regular part: V^'^ = V^^^^^^^ + V^gg^iar' such that 

limfc^o Regular = « ^nd hnifc^o A:^ ^^^guiar = ^s'l'nguiar- ^ut then we can simply redefine 
our hadronic-current as 



jj"ti ^ fj'ti + /c^K^!; , with form factors such that hm = , (4.12) 

° fc— >o 

and our initial assumption is trivially true. This proves that one can construct Eqns. 
and ([4. 3D in such a way, that V^,y is regular at A; = and f±{t') = f±{t') resulting in the 
emission graphs a) and b). The premise that limfc_j.o f± is regular can always be fulfilled by 
shifting any singular contribution into V^^- The latter must then contain a singular part 
itself and cancel the pole of f± exactly, otherwise Eqn. (J^^ cannot be fulfilled, and gauge 
invariance would be spoilt. 



The authors of [17| used another Ansatz to proof the splitting of Eqn. ( [4.7] ) with 
f±{t') = f±{t') for their studies of the radiative corrections to iCes-decays by checking 
its validity at one- loop in Chiral-perturbation theory (xPT). Given this observation, the 
authors of conclude that such a separation is a general feature of QCD. 

It remains to explicitly construct the residual IB contributions through the Ward 
identity. Using Eqn. ( [4. 7] ), the Ward identity and 

Zl{t,k) = J^f ~ (/±(t') - /±(t)) , (4.13) 
k-[pB-p) 

the missing contributions can be shown to be 

ku V^'^ = -{pB+p- kY Zf^ -{pB-p-kY Zf_ - A:^ h{t) + k, V^^ , (4.14) 

and 

ku V^'" = -{pB+P + kY Zj^-{pB-p-kY Zf_ + k^' /2 {t) + k, V^^ . (4.15) 

Using the Dirac equation, four-momentum conservation and the inverse Ward-transformation 
k'^ ^ e*^ result in the complete emission matrix elements: 

Y.M\=^e^V^^u{p^)P^v{p,)(^-^^r i^''^ + ^Pi)-£^^'iPB). + ^HPB). 



{2^ + mi) [pb]^ {mi) (pb)^ - fslu 



PB ■ k ^ pB-k 

+ (2|i + mi) Zf^ , + (mi) Zf_ , + 7^ {V^^ - Af^) j e*''{k) , 

(4.16) 
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and 



^-^ V2 \ ^Pi • p'k p'k 



r - m) Pu 7 {mi) Pv + /2 iv 

p ■ K p ■ k 

+ {2^B - mi) Zf^ , + {jni) Zf_ , + ^ {V^^ - ) e*^{k) 



respectively. Neglecting the V'^ — contributions results in an approximative descrip- 
tion of the emission process, that is exact near the soft-photon limit, where SD contributions 
do not contribute: The matrix elements Eqn. ( [4.16| ) and ( 4.17| ) then treat the mesons as 
point-like particles with structure- dependent corrections. 

4.2 IR divergencies 

Following the standard prescription of QED, we introduce a small photon mass A, that 
cancels exactly with the IR divergencies produced by the virtual diagrams. Regularizing 
these divergencies allow for the separate evaluation of the real and virtual matrix elements. 
The artificial dependence on A cancels when both are summed. Such a massive photon 
possesses a non-vanishing longitudinal polarization, what has to be taken into account. 
The three polarization vectors of a real photon emitted in the (0,0,0,l) direction are 

e^(A;) = -^(0,l,i,0), e^(A;) = -^(0, 1, -i, 0) , and e^(A;) = + A^, 0, 0, A:) . 

(4.18) 

4.3 Virtual corrections 

In this section, the short- and long-distance next-to-leading order virtual corrections are 



derived, following the splitting of Eqn. (p.4|) and the set of arguments of Sec. 3.3: the 



short-distance corrections are calculated within the Standard Model and the long-distance 



corrections modeled with the leading operator of Eqn. (2.3) and the known higher order 
tree-level operators responsible for the t dependence. The long-distance virtual diagrams 
are calculated in the on-shell renormalization scheme and at first regularized using di- 
mensional regularization. We use the standard notation of Passarino and Veltman p7[ ] to 
express the arising loop-integrals as B, and C functions, briefly summarized in App. 
The matching of the long-distance loop-integrals to the short-distance result is done via 



the Pauli-Villars prescription 1 1£ ] and the transformation of the long-distance loop-integrals 



from dimensional regularization to Pauli-Villars is further discussed in Sec. 4.4. 



4.3.1 Short-distance corrections and renormalization of Gf 



The electroweak virtual corrections to Eqn. ( p.2| ) at parton level were calculated and 
renormalized to the Standard Model by |6|, ^. The author made use of the short-distance 



(4.17) 
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expansion of the Standard Model valid at high energies, such that no leading corrections 
from the strong interaction occur. Although the latter presents a neat feature, it is beyond 
next-to-leading-order and therefore beyond the order of correction we are interested to 
examine. The relevant tree-level matrix element in [^, ^ is given by the leading operator 
ofEqn. Q, 



Ml 



with /i^ = x(g') 7m K'Z) , (4.19) 



where the hats indicate not renormalized quantities. The full electroweak corrections to 
Eqn. ( 4.19| ) involve the exchange of virtual photons, W- and Z-bosons, and Higgs scalars 
(since we are working in Feynman gauge). The electromagnetic diagrams relevant for the 
short- and long-distance photon loop-integral splitting of Eqn. ( |3.4[ ) are depicted in Fig. |^: 
Diagrams a) - c) are inter-particle photon exchange diagrams between the lepton and the 
incoming or outgoing quark-line; diagrams d) - f) depict wave- function renormalizations; 
and diagrams g) - i) are corrections due to the coupling to the VF-boson. The leading 
logarithm of the overall corrections is given by 



sd 



31n!^ + 6|Q|ln^ 



3|0|ln 



m 



w 



rrir 



+ ... 



(4.20) 



with A4q = GpA^Q, i.e. the bare leading order matrix element stripped off the bare 
Fermi coupling constant Gp; it is Q = the averaged electromagnetic charge of the 
decaying quark line, and the ellipses denote further non-leading terms. Following the 
original prescription of ^, the arising ultraviolet divergence in Eqn. (4.20) has been 
regularized by the replacement 



1 



m 



W 



A;2 - ^2 



(4.21) 



of the photon propagator, leading to a Pauli-Villars regularization, that diverges as the W- 
boson mass mw is taken to infinity. In addition a infrared regulator //q was introduced as a 
photon mass. The ultraviolet regularization in terms of the VF-boson mass in Eqn. (^ ) IS 
a priori not problematic, as long as a renormalization procedure of the bare Fermi coupling 
G-p involves a matching to a measured Fermi coupling obtained via the same regularization 
procedure, such that the explicit (but artificial) dependence in the leading logarithm of the 
VF-boson mass cancels. Such a renormalization is calculated by the same author for the 
muon decay. 



Gp = Gp 



1 + 



^In^ 

8 ml 



(4.22) 



where Gp is the measured Fermi decay constant. This leads to the desired short-distance 
virtual matrix element 



sd 



^--(l + 2|Q|)ln- 



+ ... , 



(4.23) 
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where the eUipses denote further non- leading terms, and the leading matrix element A4q 
contains the renormalized Fermi coupling Gp . The necessity of the involved virtual photons 
to resolve the involved mesons of the decay usually leads to the choice of a hadronic scale 
for /xq: a choice often used for semileptonic S-meson decays is given by the 6-quark mass 
in the MS scheme: /ig = i^b = ^h^J^bl- The logarithm in Eqn. ( f4.23| ) represents then the 
short-distance corrections due to virtual particle exchange ranging from energies between 
mf, and mz- In order to obtain the desired separation of long- and short-distance scales 
in our model, we choose a lower matching scale and fix //q = "^n^x for all decay modes. 
The logarithm in Eqn. ( 4.23| ) represents then the short-distance corrections due to virtual 



particle exchange ranging from mx to mz- 

4.3.2 Next-to-next-to-leading order corrections at ©(aQsGp) 

QCD loop corrections change the form of the effective Hamiltonian of the matrix element 
Eqn. ( p^ to 



2 

I 2^ 



[i'vl^Pi.UK^ (4-24) 



with C 



,as{fJ.'^) = 1 + O (as(/i^)). At one-loop level, assuming three colors and six 



flavors, the running of the strong coupling constant is given by as{fM^) = 47r/ ( /3o In -r^ — 
with /3o = 7. Then the correction due to strong interaction effects is approximatively |^ 

^o,sd = ^lO|/a.K> (4.25) 



with 

dfc2 — a,(A;2)) . (4.26) 



Mo 



With m\Y = 80.34 GeV, Aqcd = 0.2 GeV and /xq = w-b = 4.2 GeV this leads to a correction 
of —0.03% to the transition probability predicted by Eqn. ( [4.23| ). Therefore for /xq = rrib 
strong correction effects, formally at next-to-next-to-leading order and 0{a Gp) turn 
out to be completely negligible. The estimated corrections at 0{Gp aas) analogous to 
Eqn. (|]2|) for fiQ foi B Dlu (7) and B D^lv (7) are summarized in Table |^. Note 
that one is free to choose a different matching scale ^0 for such corrections. Therefore non- 
perturbative corrections at next-to-next-to-leading order for S — >• vr / (7) with /xq = to-tt 
are not troublesome . 

4.3.3 Renormalization of masses and electromagnetic coupling 

The corresponding long-distance correction to the short-distance result of Eqn. ( [4.23 ) is 



calculated in the following sections. We will work in the on-shell renormalization scheme, 
where the poles of the propagators are given by the phenomenological masses. In addition, 
we assume, that the electron charge is renormalized to its phenomenological measured value 
in a similar manner as the Fermi coupling constant. We will not work out the details of 
these renormalization steps, but only consider the diagrams that contribute to the virtual 
corrections of the B ^ X lu decay amplitude, as shown in Fig. (0). 
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Mo as(Mo) AO{aasGF) 



B+ ^ D^lv (7) m£,o 0.40 -0.04% 
D+lu (7) m£,+ 0.40 -0.04% 
B+^Dl^lv{-i) m^,o 0.36 -0.04% 



B'^ D*+ 1 u {-f) m^.+ 0.36 -0.04% 



^0 

^0 

Table 1: The choice for the infrared regulator /iq for 5+ — > X'~' I v (7) and — > ^ (7) decays 
with X = D OY Dq are hsted. The value for asf/io) was extracted with a^f/i^) = An/ (/3o In -r# — ) 

where mw = 80.34 GeV and Aqcd = 0.2 GeV. AC'(aas Gf) is the correction in percent to the 
pure short-distance electroweak correction: all corrections due to 0[oiOis Gp) effects are negligible 
small. 



4.3.4 Structure-dependent corrections to virtual diagrams 

For the long-distance loop calculations the meson will be treated as entirely point-like, 
neglecting all structure-dependent related corrections beyond the gauge fixing terms needed 
to fulfill the Ward- identity. Their coupling can be isolated in Eqns. ( 4.16| ) and ( [4.17 ) by 
setting Zf^ = and their contributions need to be included to obtain a consistent result 
with respect to gauge invariance. As long as ^ mx this assumption is justified, and the 
long-distnace loop integrals only receive additional corrections from higher order resonances 
through B*/**/-X'^ and X*/**/---X'^ coupling, which we neglect. 
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a) 4 b) c) 




d) e) /) 




Figure 4: The electromagnetic short-distance one-loop photon diagrams for 6+ — > x^^ Iv and 
^ x~lu are shown: a) - c) inter-particle photon exchange; d) - f) wave-function corrections; 
and g) - i) coupling from external legs to the PF-boson. 



a) b) 




f) 9) h) 




Figure 5: The long-distance one-loop diagrams for B+ I u and — >■ X~ I v are shown: 

a) and h) inter-quark photon exchange; c) - e) wave- function corrections; and f) - h) structure- 
dependent corrections. 
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4.3.5 Long-distance wave function renormalization 

The fermionic long-distance wave function correction, Fig. ^c), is given by 



1 d^i^i 



2 d^i 



1 a 
24^ 



1 + 2Bi{pf- mlX^)+ 4mfBi{pf;ml X') + 8mfBo{pf;ml A 



2.^2 \2\ 



(4.27) 



where S(|^;) is the fermionic self-energy correction. The long-distance wave function cor- 
rections of the scalar legs, Fig. ^ d) and e), are given by 



1 dS _ 1 a 
2d^ ~ 24^ 



2^0 (p^ ; A^) + {4m^ - \^)Bq{p'^- m^,X^' 



(4.28) 



and 



1 dE la 



2 dpj^ 2 47r 



2Bo{pl;ml,X^) + (4m| - A^) (p| ; m| , A^ ^ 



(4.29) 



respectively, with S the scalar self-energy correction. The total correction due to long- 
distance wave function renormalization is then given by 



Ml 



d^i) ^ dE 
d|5; dp"^ 



' 



(4.30) 



or 



d^iM I dE 
d^i dp| 



Ml 



(4.31) 



The corresponding short-distance diagrams are given by Fig. - f). 



4.3.6 Long-distance inter-particle exchange 

The long-distance inter-particle exchanges. Figs, ^a) and b), result in the following matrix 
element 



Ml = (+e2)i=GFKb [u{pu)PR7^.r''v{pl)] 



(4.32) 
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with the correction tensor as either 
i 



(47r)2 
i 



(47r)2 

i 



+ 
+ 
+ 



(47r)2 
i 

i 

(4^L 



(/3(i) Pf) (4pB • + (if^) [4pB ■ PI + 2mf - 2^B^i))]Ci{pf , s,pl; X^mf, m%) 
(/3(i)p^) (4PB • Pi) + {H^) {2pB ■ PI + 2ml)]c2{pl s,ph A^, mf, m|) 



Mt)p';)UpB-pi + 2mf -2^B^i) Cn{pf,s,pl;X^mlml) 



Mt)PB)['^PB ■Pi + 2m%] C22{pf,s,PB'^X^,mf,m%) 



+ 
+ 



(47r)2 



(47r)2 
i 

(4^ 



Mt)pf)(2pB-pi + 2ml 



+ - f3{t)p''B)UpB-pi + 2mf -2^B^i) Ci2{pt,s,pi;X',mf,m%) 



-Mt))[2pf + 2^''^B)\Coo{pt,s,p'B;X',mf,mj,) 
H^' + h{t)pf\Bo{s-mlml) + ^\f^it){pi - pBY]B^{s;mlml 



(4.33) 



with s = ipi — pb) , or 



+ 



+ 



(47r)2 
i 

(4^ 

i 

(4^ 
i 

(4^ 

i 

(4^ 
(4^ 



- 4(p • Pi)) J Coipf, s,p^- X\ml w?) 
f2{t)p';){-4p-pi) + {H>')(4p-pi- 2mf - 2ii>i))'\ Ci {pi s,p^; X^ ,ml w? 
{f2{t) P^) (-4p • pi) + (2m2 - 2p ■ pi)]C2{pl s,p^; X\ml n?) 
f2{t) ) (Ap • pi - 2mf - 201^ ] Cu {pI s, p^; Xlmf,m^) 
f2{t) p^) (2m2 - 2p ■ pi)] C22{pf, s,p^;X'',mf,m^) 



(^f2{t)pf) (2m^-2p-pi 



+ ( f2{t)p^ )Up-pi-2mf- 2i,i>i ) Ci2{pi, s,p^;Xlmlm^) 



(47r)2 L 
i 



/2(t) j [-2p^ + 2r^)\Coo{pf, s,pl X\mt m') 
^^ + /2(t)pf]So(s;m2,m2) + ^[/2(t)(pi+p)^]i?i(s;m2,m2), 



(4.34) 



with s = {pi+ p) , respectively. The corresponding short-distance diagrams are given by 
Fig. I a) - c). 
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4.3.7 Structure-dependent contributions 

The long-distance one-loop diagram between the remaining structure-dependent contribu- 
tions in the IB coupling and the external legs, Figs. ^ / - /i, are given by 

o.Gf 



Ml 



V2 



(4.35) 



with the correction matrix A, which is either 



A 



(47r)2 

i 



1 r 



(47r)2 mf 
i 



(47r)2 
i 1 



(47r)2 2m|j 



[2i>B h{t)\Bo{ml-X\ml) 
[- i>B h{t)\Ao{ml) 



(47r)2 



or 



A 



(47r)2 
i 



-2i>if2{t)\ Bo{ml,X^mf) + 



(4.36) 



(4vr) 



+ 
+ 



1 r 



(47r)2 
i 



(47r)2 2m2 



[2i,f2{t)\B^{n?-\\m^) 



(47r)2 



i>if2{t) 



(4.37) 



The corresponding short-distance diagrams are given by Fig. ^g) - i). 
4.4 Pauli-Villars regularization 

It remains to transform the long-distance loop-expressions from dimensional regularization 
into the Pauli-Villars regularization scheme. This can be done by a simple set of replace- 
ment rules and counter terms. The Pauli-Villars prescription introduces a massive 
photon field with opposite norm into the theory via 

1 



(4.38) 



Taking the limit of /Uq — )• oo decouples the field from the rest of the theory and gives raise 
to a UV divergency growing as ~ Ih/Uq, which is exactly the desired opposite behavior of 
the short-distance result. The massive field suppresses the exchange of virtual particles 
with energies larger than ^o- 

4.4.1 Wave function counter terms in Pauli-Villars 

Eqn. ( [4.27 ) has to be modified via 

Bo{mf,mf, A^) Bo{mf;mf, A^) - Bo{mf;mf,fil) , 

Biimlmf, A2) ^ Hmlml A^) - Bi{mf,mf, fi^) , (4.39) 

with i = 0, 1. Adding the counter-term 

dS(|^i)pV terms 



d|5; 



a 
4tt 



(4.40) 
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results in the desired transformation. Similarly, Eqn. ( 4.28| ) replaces 



Bo{m^,m'^, A^) BQ{m^,m^, A^) - Bq{w? ,m? , ixq) , 



with a counter-term 



dp' 



PV terms _ 



(4.41) 



(4.42) 



The rules for Eqn. (|]2|) can be obtained by replacing m? — )■ m\ in Eqns. ( [4.41] ) and 



(14421). 



4.4.2 Inter-particle exchange diagrams in Pauli-Villars 

Eqn. ( [4.33 ) has to be modified via 



Ci{pf,s,p%;X'^,mf,m%) Ci{pf , s,p%; ,mf ,m%) - Ci{pf , s,p%; fil,mf ,m%) , 
Cij{pf,s,p%;X'^,mf,m%) djipf , s,p%; ,mf ,m%) - djipf , s,p%; ^xl,mf ,m%) , 

(4.43) 



Bi{s;mf,m%) 0, 

with i, J = 0, 1, 2. Adding the counter terms 

i 



PV terms 



(4vr)^L 



+ - 



+ - 



(47r) 

i 



47r 



H'' ^oj Co{pf, s,p%; ^0. "if , m|) 
fsit) {piY l^l Ci{p^,s,p%;nl,m^,m%) 
h{t) {pb)'' fJ'l C2{p'i,s,pl;nl,mf,m%). 



(4.44) 



results in the desired transformation. For Eqn. ( [4.34[ ) one modifies 

r 2 2 \2 2 2\ > / 2 2 \2 2 2\ r~< / 2 2 2 2 2\ 
Ci{pi,s,p ;A ,mi,m ) Ci{pi,s,p ;A ,m;,m ) -Ci{pi,s,p ■,ijQ,mi,m ), 

r< I 2 2\2 2 2\>/~(/2 2\2 2 2\ / 2 22 2 2\ 

Cij{Pi,s,p ;A ,m;,m ) ^ Cij{pi,s,p ; X ,mi,m ) -Cij{pi,s,p ;/io,"i,,m ), 



0, 



(4.45) 



with i, j = 0, 1, 2 and adds the counter terms 

H>' fil Co{pf,s,p'^;fil,mf,m%) 



PV terms 



+- 



(4vr)' 

i 

(4^ 



/2(t) (w)^/^o Ci{P^,s,p^;^j.l,mf,m^) 
f2{t) {pY fil C2ipf,s,p^;fil,m^,m^) . 



(4.46) 
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4.4.3 Structure-dependent contributions in Pauli-Villars 



Eqn. ( 4.36 ) is modified according to 

Bi{mf;X'^,mf) Bi{mf;X^,mf)-Bi{m^;^j.l,mf), 

with i = 0, 1. Adding the counter-term 



Apv 



terms 



(47r)2 



(4.47) 



(4.48) 



results in the desired transformation. Eqn. ( [4. 371) has to be modified according to 

Bi{mf; X^,mf) Bi{mf; \^,mf) - Bi{mf; fJ-l,mf) , 
Bi{rn^; A^, m^) — )• Bi{m'^; )? ,m^) — Bi{m'^; fi^jin^) , 



with i = 0, 1, and a counter term of 

Apv terms 



(47r)2 



(4.49) 



(4.50) 



4.5 Next-to-leading order difTerential and total decay rate 

Summing all virtual long-distance corrections for B^ — t- I v or B^ — t- I u, given by 
Eqns. ( 4.32D , and either ( 4.30 ) or ( [4.31 ), and ( 4.35 ), and the short-distance result Eqn. 
(4.23), yields the approximative virtual next-to- leading order differential decay rate. It is 
given by 



dro + dr^ 



1 



64 Ti-^mB 



Ml\ +2^MlMl + 2\Ml 



0|2 



r3 a 



^-(l + 2|g|)ln^ 
4 7r' ^0 




Similarly the real corrections with Eqn. ( 4. IE ) or ( 4.17] ) result in the approximative real 
next-to-leading order differential decay rate: 



drf 



{2Tiy^EEiE,Ek 



5*^^) {niB -p-pi-Pu-k)\^^Ml\ d?p d^pi d^p,y d^k , 



(4.52) 



with Ek = and E,y = P^- Integrating Eqns. ( 4.51 ) and ( 4.52| ) yields the approxi- 
mative next-to-leading order total decay rate. Comparing with the total tree-level decay 
rate, the integral over phasespace of Eqn. ( |2.7D , yields the correction factor (5totai due to 
next-to- leading order effects: it contains all short- and long-distance next-to-leading order 
corrections. It is 



1 + 5, 



total 



J- 



l + 5sd + 5ld 



(4.53) 



where (5total = ^^sd + ^id, with the short-distance contribution 5sd = — In ^ from Eqn. 



( 4.23 ), and 5id denotes the long-distance corrections. 
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5. Numerical evaluation 



In order to evaluate Eqns. ( 4.51| ) and (4.52), we developed the Monte Carlo generator 



BLDR {B-meson leading order radiative corrections), which is derived from the code of KLOR 
{Kaon leading order radiative corrections - written by Troy Andre ||l^). It incorporates all 
form-factor models of App. ^ and the next-to- leading order matrix elements derived in the 
previous section. The loop-integrals are evaluated numerically using the software package 



Looptools |28] and the regularization procedure can be changed between dimensional 
regularization, Pauli-Villars and an approximative Euclidian cut-off method. The phase- 
space integration occurs in two steps: the total tree- level and next-to- leading order rates 
are calculated to determine 5\d and the IR cut-off dependent radiative fraction of: 

pO I T-l 

^ + 



(5.1) 



This is done using the VEGAS algorithm described in [^] . Subsequently the radiative frac- 
tion Eqn. ( |5.1| ) is used to mix events generated from the real and virtual next-to-leading 
order matrix elements using standard Monte Carlo techniques. In order to optimize the 
event generation procedure, moderately pre-tuned probability density functions of the in- 
tegration variables in question are used. In addition to the implemented next-to- leading 
order calculations, the BLOR Monte Carlo generator can be interfaced with the approxima- 



tive next-to-leading order package PHOTOS The code of BLOR is publicly accessible 



and can be found at |30|]. The particle masses and couplings implemented therein are 
summarized in Table 0. 



Particle 


Mass 


mx(45) 


10.5794 GeV/c^ 


rT(4S) 


20.5 MeV/c^ 


mB+ 


5.27913 GeV/c2 


m^o 


5.27950 GeV/c^ 




1.86950 GeV/c2 




1.86484 GeV/c2 




0.13957 GeV/c2 




0.13498 GeV/c2 


^0 


2.40300 GeV/c2 




2.35200 GeV/c^ 



Parameter 



Value 



A 

a 



10-7 GeV/c2 
0.00729735039 
1.16637 • 10-5 (^^GeV- 



Table 2: Particle masses and couplings used in the simulation. All values are taken from fell 
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6. Predictions 



In the following, we quote the value of (5totai calculated from the next-to-leading order 
model and try to assess its uncertainty due to the approximative character of the matching 
procedure and missing real and virtual SD contributions for B ^ Dlu (7), B ^ ttIv (7), 
and B — >• Dq I v (7). Then, we present the main result of this study: normalized kinematical 
distributions for the final state lepton, meson, and photon. For a comparison of our findings 



with the next-to-leading order algorithms of PHOTOS and PHOTONS++ [32|, we refer the reader 



to |3| 



6.1 Uncertainties of Jtotai 

Uncertainties to the next-to-leading order corrections come from five sources: the numeri- 
cal evaluation ((Tnumericai), the matching procedure between long- and short-distance result 
(^matching)! the limitations of the effective theory describing the complete phase-space ade- 
quately due to neglecting the V^,? — A^]^ contributions ((Tsd)i next-to-next-to-leading order 
electromagnetic effects ((Tnnio)i and last but not least model dependent uncertainties. The 
total uncertainty by the first four sources is given by 

"^total — "^numerical "r "^nnlo + "^matching + "^SD ) \^-^) 

It is not possible to assess the model dependent uncertainties due to the negligence of the 
unknown higher-order operators of Eqn. (^]^), therefore the stated uncertainties should 
be interpreted as uncertainties that can be assessed within the model itself. Without full 
knowledge of the full effective theory, no deeper and more satisfying assessment can be 
made. In the following the four sources of uncertainties that can be assessed are discussed. 

6.1.1 Numerical uncertainty of the integration 

The total integration error of <5totai is given by 

m ^ 1 



"^numerical — ^ ^ (j? ) ' (^■^) 



where fj, is the error of the zth VEGAS evaluation, and m the number of total evaluation 
steps. 

6.1.2 Corrections due to next-to-next-to-leading order electromagnetic correc- 
tions 

The next-to-next-to-leading order (nnlo) electromagnetic corrections are estimated as 

o-nnio = a ( (^sd + (^id ) • (6.3) 
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6.1.3 Uncertainty due to the matching procedure 

The uncertainty due to the matching procedure is studied by varying the matching scale 
yUo within reasonable bounds. By doing so, one varies with what contributions of the 
effective Hamiltonian Eqn. the 'intermediate region' with respect to the short- and 

long-distance regime are modeled. This unveils a shift in 5totai = <^totai(/^o) and we choose 
to vary the matching scale within ^/xq and 2/xo. The resulting estimator for the matching 
uncertainty is then calculated as 

^matching = ^ I '^total (A'o) " ^totaK^^o)! + ^ I '^total (/^o) " '5total(2Mo) | • (6.4) 

Note that for a matching scale choosen at 2/xo, the phenomenological model for the long- 
distance contributions is beyond its expected validity; whereas fixing the matching scale 
at ^//o short- and long-distance corrections should represent a reasonable description. The 
dependence of the matching scale fiQ on the overall correction ^totai are summarized in 
Tables 0, and |l3. 



6.1.4 Corrections due to real SD contributions for B ^ D I u j decays 

For B — )• Dlu^ decays, the dominant corrections from hadronic emissions lie in the soft- 
photon part of phase space and are due to IB contributions. The additional SD contribu- 
tions in this region of phase-space are well approximated by emissions from the two lowest 
excited states, B* and D*, respectively. These intermediate emissions are depicted in Fig. 
(|6D and their contributions to the SD vector and axial-vector coupling are 



[p + kj — rrij-ft [P — K) — m^, 

(6.5) 



_^SD _ i{D\Jr\D*){D*\y,-WB) ^ i{D\Vf.-A^\B*){B*\jr\B 



where the ellipses denote contributions from higher excitations and further correction terms. 
The emission from the L'*-line is expected to be dominant over i?*-pole contribution, since 
the latter can only occur off-shell and the large mass of the B* excitation suppresses 
emissions considerably. The intermediate D* is kinematically allowed to be on-shell in 
certain regions of phase-space. As proposed by Q in these on-shell regions one can use the 
lattice D*Dj coupling of |Q to estimate the first term of Eqn. ( |6.5D as 



2p-k + Am|j, 

- {pB+P + k)^{pB-p-k)^a+{t') 

- {PB-p-k)^{pB+p + k)^a_{t')) , (6.6) 



where gD*D-y denotes the corresponding effective D*Dj coupling, Am|,* = m'j-, — m^* + 
iniD* ^D* with To* the D* width, and /, g, and a± the on-shell B — t- D* form factors, 
further discussed in App. Assuming the dominant contributions occur near the on- 
shell value of the propagator, the additional SD contributions to the hadronic B ^ Dlv^ 
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emissions can be estimated by including Eqn. (6^) into the emission matrix element. The 
resulting shifts in Jtotai are summarized in Table ^, where Jtotai + SD approx. denotes the 
total predicted correction gained by including Eqn. ( [6.61) . The difference to (5totai is taken 
as uncertainty ctsd due to real SD emissions: 

<7SD = I'^total — (^total + SD approx. I • (6.7) 

Contributions from higher scalar and tensor charmed resonances at 0(Gfq;) to real B — )• 
Dlvy decays are forbidden by parity, angular momentum and spin conservation. Higher 
vector charmed resonances, however, can contribute at 0{Gpa). 




Figure 6: The SD one-particle contributions to the real emission are shown: a) and b) emission 
from a B*Bj and D*Dj vertex, respectively. 

l^total (^total + SD approx. C^SD 

^D+lv (7) 0.02223(6) 0.02225(7) 0.00002 
B+^D°lu{j) 0.01463(5) 0.01627(6) 0.00158 

Table 3: The effect of including the approximative SD emissions on S\d ioY B ^ D ev (7) is shown. 
The difference between the charged and uncharged charmed final states originates from the effective 
coupHng go'D-i- it is go'+D+j — —0-1(7) and gD">D"-y = 2.7(1.2) from [ p4| . 



6.1.5 Corrections due to real SD contributions for B ^ it 1 1/ ^ 

For B ^ irluj the real IB result can be compared with the findings of [^]: The authors 
calculated the real IB and SD corrections in the framework of SCET and in the phase-space 
region of soft pions and hard photons: E-,^ < 0.5 GeV, k > 1.0 GeV, and 0^^ > 5". In 
this region a significant influence of SD contributions can be expected, since the first light 
intermediate resonance receive on-shell contributions at fc ~ 0.8 GeV. Since the splitting 
of in IB and SD contributions differs from our construction by not incorporating 
corrections of 0{k) and beyond into their IB term, the comparison also offers interesting 
insight into the predictive power of including such terms. In order to distinguish this we 
write IB' in the following when referring to their choice. The calculation of |25| uses the 
leading order heavy hadron chiral perturbation theory (HH^pt) form factor predictions: 

ME.) oc ±^ (6.8) 
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with E'tt = , and A = niB* —niB = 50 MeV. The radiative form factors in our formulation 
of the IB emissions are thus given by f±{E.j^ + Ek) for a charged final-state pion, whereas 
in the IB' choice, Eqn. (^^) is used as presented. The partial and differential branching 
fractions from both calculations are compared in Table Q and Fig. |^: the pure IB based 
prediction based on the model dependent choice in this work, reproduces the photon energy 
spectrum well. The authors of |25| also quote the results of the partial branching fraction 
gained from the IB' terms alone, and from including the leading terms from Low's theorem 
{k~^ and k^), denoted as IB" in the following. The partial branching fractions for the IB' 
and IB" calculations are stated in Table |5|, and differ significantly from the IB'+SD and 
the IB predictions. 



k energy: — > tt^ ei'j 

T 0.09 
> 

O 0.08 

0.07 
"b°0.06 

0.05 

0.04 

0.03 

0.02 

0.01 


1 1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 

k [GeV] 

Figure 7: The differential branching fractions of the B'^ 7r+ e 7 decay is shown: blue shows 
the result of ||2^, black is the prediction based on the IB contributions with HHxpt form factors, 
and grey is the prediction based on the IB contributions with Ball-Zwicky form factors. 






Cuts 


IB' + SD HHxpt [H 


IB HHxpt IB Ball-Zwicky 




> 1.0 GeV 






EjT 


< 0.5 GeV 

> 5° 


1.2 1^;^ X 10-6 


1.16(1) X IQ-^ 0.35(1) X 10-6 



Table 4: The partial branching fractions BF{B'^ 7r+ ev^)cut for the IB'+SD and IB calculation 
are stated. The errors in the parentheses are numerical. 



Therefore including terms of O^k"^) and beyond into the IB contribution via Eqn. 
(4.17) reproduces the main features of the photon energy spectra and provides a partial 
branching fraction close to the result of the full IB' + SD calculation. 

By varying the HHxpt IB result of the partial branching fraction to the extremal values 
within the error band of the result from |25|, one obtains the estimator for the uncertainty 
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Cuts 










IB' HHxpt 12§ 


IB" HHxpt ^ 


E-y 


> 1.0 GeV 








< 0.5 GeV 


2.8 X 10"^ 


2.4 X 10^6 




0e7 > 5" 







Table 5: The partial branching fractions BJ-{B^ 7r+ ej^7)cut for the IB' and IB" based predic- 
tions. 

(TsD due to SD effects in ^ I v ^ in this region of phase-space, as 

0"SD = l^total IB - ^total IB + Sd| = 0.0048 , (6.9) 

We interprete this as an estimator for the uncertainty for the entire phase space. Comparing 
the HHxpt prediction of (5totai and the partial branching fraction with the result obtained 
with the Ball-Zwicky form factors (see App. reveals a deviation due to their different 
parametrization and normalization. The predictions for (5totai sue listed in Table |6|, and the 
differential and partial branching fractions are depicted in Fig. and listed in Table |^. 
The deviation in ^totai between HHxpt and Ball-Zwicky is small with respect to (Tsd- This 
suggests that Eqn. (|6.9| ) can be used as a valid estimator for the missing SD contributions 
for the ^totai prediction based on the Ball-Zwicky form factors. 

B° -^TT+eu (7) HHxpt Ball-Zwicky 
~ I-Fi^sd + i^id 0.0562(1) 0.0555(1) 

Table 6: The prediction for l-t-5sd+<5id for HHxpt and Ball-Zwicky form factors for ^ 7r+ e (7) 
are shown. The uncertainties in the parentheses are numerical. 

For the radiative — t- lu^ decay no calculations exists that estimates the real SD 
contributions, therefore no uncertainty due to such corrections can be assigned to our IB 
based prediction. 

6.1.6 Corrections due to real SD contributions from B ^ DqIvj 

For B ^ DqIu^ decays significant corrections due to SD contributions can be expected: 
the short lifetime of the Dg-meson results in a large decay width, which causes an overlap 
with the broad charged Dl~^- and the narrow uncharged -Djj' '^-states. We are not aware 
of any lattice results, which examine the coupling between these two charmed meson res- 
onances, the DQ-meson and the photon, or of any work that could help estimating this 
uncertainty. 

6.2 Predictions for (5totai 

The predictions for 1 -|- (5sd + (^w and + 5sd + ^id for B — Dlv{'^) are listed in 

Table |^. The assigned uncertainty interval estimates possible further contributions due 
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to real SD terms, matching uncertainties, numerical uncertainties, and next-to-next-to- 



leading order corrections - within the limitation addressed in Sec. 6.1. The mass of the 
D— meson allows a reasonable separation of short- and long-distance domains. Table ^ 
summarizes the predictions of 1 + 4d + ^id and 1 / \/l + 6sd + <Jid for B — )■ tt lu (j). The 
assigned uncertainty interval estimates the matching uncertainties, numerical uncertain- 
ties, and next-to-next-to- leading order corrections. In addition for B^ — )• vr^ / v (7) the 
real SD corrections were estimated. Table ^ eventually summarizes the predictions of 
1 + 5sd + 5id and l/\/rT4dT^ for B ^ D*oli^ (7). The assigned uncertainty interval 
estimates the matching uncertainties, numerical uncertainties, and next-to-next-to-leading 
order corrections. In addition the mass of the Dq— meson allows a reasonable separation 
of the short- and long-distance separation. But current knowledge does not allow to make 
any assessment for real SD corrections. 

The inverse of y/1 + 6sd + ^id is the relevant correction to the CKM matrix elements 
l^cbl and l^ubl gained from measured decay rates. In summary, the source of the domi- 
nant uncertainty in our presented approach originate the matching procedure. In addition 
unknown uncertainties due to the neglegence of unknown SD contributions and model 
dependences arise. 

The predicted corrections break as expected the isospin symmetry of the semileptonic 
decay process, hence the dominant real and virtual corrections scale with the mass of the 
charged meson in the decay. 

1 + 5sd + Sid 1/Vi + Ssd + 5id 





D- 


" (7) 


1.0222(1 ± 2± 17± 1) 


0.9891(1 ± 1 ±4± 1) 


B^- 






1.0222(1 ± 2± 17± 1) 


0.9891(1 ± 1 ±4± 1) 


B+ 




" (7) 


1.0146(1 ±1 ±39 ±16) 


0.9928(1 ±1± 10 ±4) 


B+ - 






1.0147(1 ± 1±39± 16) 


0.9927(1 ±1± 10 ±4) 



Table 7: Predictions for (5totai ioT B ^ Dlv (7) are listed. The uncertainties in the parentheses 
are numerical, next-to-next-to-leading order, matching, and due to neglected real SD contributions, 
respectively. 

1 ± <5sd + 5id 1/Vl + 4d ± 



B^- 


vr e+ z^e (7) 


1.0555(1 ±4± 148±48) 


0.9734(1 ±1 ±33 ±10) 


B^- 


> TT" /i+ (7) 


1.0545(1 ±4± 136±48) 


0.9738(1 ±1 ±31 ±11) 


B+ 


-7- vr'^ e"*" V(, (7) 


1.0411(1 ±3± 100) 


0.9801(1 ± 1 ± 23) 


B+ - 


vfO 11+ (7) 


1.0401(1 ±3± 89) 


0.9805(1 ± 1 ± 21) 



Table 8: Predictions for (Jtotai for B ^ ttIv (7) are listed. The uncertainties in the parentheses 
are numerical, next-to-next-to-leading order, matching, and for _B" tt^ Iv (7) due to real SD 
contributions, respectively. 
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1 + '^sd + Sid i/yi + '^.d + 

S° ^ " e+ i/e (7) 1.0224(1 ± 2 ± 10) 0.9890(1 ± 1 ± 2) 

^ D*- fi+ i^^i-f) 1.0226(1 ± 2 ± 10) 0.9889(1 ± 1 ± 2) 

B+ ^ D*Q° e+ Ve (7) 1.0142(1 ± 1 ± 35) 0.9930(1 ± 1 ± 8) 

B+ D*o fi+ (7) 1.0144(1 ifc 1 ifc 35) 0.9929(1 ifc 1 ifc 8) 

Table 9: Predictions for (5totai for B Dq I v (7) are listed. The uncertainties in the parentheses 
are numerical, next-to-next-to-leading order, and matching, respectively. 



6.2.1 Corrections to |V"cb| 

The lepton mass for / = e or only has a negligible effect on the value of the total decay 
rate. One can therefore assume that 

r[!(B+^X°ez^) =r[!(5+^XVz^) and r[!(B°^X+ez^) =r[!(50 ^X+^i^), 

(6.10) 

and average the results of 1 -|- b^^ + S\A and \/ 1 -|- (^gd + S\A iii Table ^ from both lepton 
channels to decrease the numerical uncertainty of the integration. The SD and next-to- 
next-to- leading order uncertainties then are correlated by 100%. This yields an averaged 
correction of 2.22(17)% to the total — t- D^lv decay rate in comparison to the tree- 
level decay rate. The increased total decay rate translates into a decrease of the extracted 
value of 1 14b I from measurements by —1.09(4)% in comparison to extractions based on the 
tree-level calculation. Similarly, the B^ — )• Z)*^ I v total decay rate increases by 1.46(43)%, 
translating into a decrease of —0.72(10)% of |Fcb|. 

Analyses often study S-mesons from the decay chain T(4S) B B with B either B^ 
or B^ . In order to reduce the statistical uncertainty, and assuming isospin, the B^ — J- Tf' I v 
and B^ — J- / v decay rates can be combined. At tree-level, this combined total decay 
rate is given by 

)^1{B ^Dlv) = f^_ r[](i?+ ^ DO / 1/) + /oo r[](i?o ^D^lv), (6.11) 

with 

/+_ = i3-F(T(45) ^ 5+ B-) and /oo = BT(T[^S) B^ B^) , (6.12) 

with the branching fractions of the T(45') resonance either decaying into a B^ B^ or B^ B^ 
pair. Including the isospin conserving parton-level correction in Eqn. ( [4. 231) yields 

(i + 4d)r[J(i3^i?/i^) « l^cbl', (6.13) 

with 5sd = 1.430%, resulting in a correction of —0.707% to the extracted value of |V"cb|. 
Including the isospin breaking contributions from the long-distance corrections of Table 
result in 

(1 + 4d + <5i^) /+- r[J(i?+ ^ 50 / z.) + (1 + 4d + /oo rl](i?o ^ / 1.) oc iFcbP , 

(6.14) 
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where 6^ and (5{^ either denote the long-distance next-to-leading order correction of 
D^l V or —7- I V decays, respectively. The total correction is 



(1 + <5total) r[l(i? ^Dlv)^ l^ebl' , (6.15) 



with 



_ Vl{B^^D^lv) Tl{B^^DHv) 



5.16) 



Assuming isospin symmetry at tree-level and averaging over both lepton modes, and 
with = 1.065(26) from [31| yields an overall correction of (^totai = 1-86(29)% due 



to short- and long-distance next-to-leading order effects^. This translates in a reduction of 
—0.91(7)% to the extracted value of \Vch\- 



1 + <^sd + 6id 1/Vl + 4d + 5id 



B^^D-lv 1.0222(17) 0.9891(4) 

B+^D^lv 1.0146(43) 0.9928(10) 

B-^Dlu 1.0186(29) 0.9909(7) 

1 + ^^.d + ^Id l/\/l + ^sd + '^Id 

B^ -^TT-lu 1.0550(150) 0.9736(34) 



Table 10: Averaged integration results for 1 + i5sd + ^id and V 1 + <5sd + iox B ^ Dlv and 
— >■ 7r+ / v. The uncertainties in the parentheses are the sum of numerical, next-to-next-to-leading 
order, matching, and due to missing real SD contributions. The result for B'^ ^ tt^ Iv should be 
used with care. 



6.2.2 Corrections to iKibl 

Averaging the leptonic results yields for the B^ — )■ vr"*" / u total decay rate an increase by 
5.50(1.50)%, translating into a decrease of —2.64(34)% of extracted values of \V^\j\ from 
measurements, as summarized in Table |l^ 

6.3 Predictions for next-to-leading-order differential rates 

In the following, the normalized differential decay rates for the tree-level and next-to- 
leading order decays for the lepton, meson, and photon are shown: 



*The uncertainty is propagated assuming a correlation of 100% between the SD and next-to-next-to- 
leading order uncertainties of each channel and a correlation of 100% between the isospin doublet. Further 
the uncertainties of /-| and /oo are correlated by —100%. 
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with r = Fq or r = Fq + Fq + . All decay rates are presented in the experimentally 
important rest frame of the T(4S') B B decay. 

Fig. 1^ depict the results for B^ — )• l)+ e v (7) and B^ — )• jjL v (7) decays. Fig. ^ 
depict correspondingly the B^ — )• eu (7) and B^ — )• fii^ (j) decays. The very mas- 
sive L'-meson only receives small corrections due to radiative effects, whereas the electron 
three-momentum distributions are notably shifted towards smaller three-momentum. The 
correction on the muon three-momentum distribution are smaller, due to its larger mass 
in comparison to the electron. 

Fig. |l^ depict the predictions for the B'^ — )■ tt'^ ei' (7) and B^ —?• tt^ (7) decay 
modes. The charged pion final state radiates considerably and receives corrections in its 
shape. This is caused by the dependence on the radiative four-momentum transfer squared 
i' = {pb — P — k^'^ the IB terms. The electron three- momentum distributions are shifted 
towards lower momenta, and the muon three- momentum distributions receive small correc- 
tions. Fig. 11 show the B^ — t- tt'^ 62^(7) and B~^ — ?■ nvlj) predictions. The uncharged 



pion receives negligible corrections from the radiative four-momentum transfer squared 
i' = {pb — p — due to the large mass of the 5-meson, that suppresses considerably real 
IB next-to- leading order corrections. 

Figs. present the predictions for the B^ — J- Dq ~^ ev (7) and B^ — ?• Dq ~^ fiu (7) 
differential decay rates. The three-momentum of the D^* , similarly as the case, only 
receives small corrections, whereas the electron three-momentum distribution is shifted to 



lower momenta, and the muon distribution receives small corrections. Similarly, Fig. |13 
depict the predictions for the B^ — )• Dq^ ev {'y) and B^ — )• Z)q 1/ (7) differential decay 
rates. 



7. Conclusion 

The next-to-leading order corrections for exclusive semileptonic i?-meson decays into pseu- 
doscalar and scalar final states were calculated in our model approach for B ^ Dlv (7) , 
B — )• 7r/z^(7), and B — J- DQlv{'y) decays. We predict the previously unknown total en- 
hancement to the total decay rate ^totai for B — )• D lu {'y): averaged over lepton species 
the B ^ Dlv (7) total decay rate receives a global correction of 1.86(29)% in comparison 
to the tree-level prediction. This results in the reduction of the extracted value of \Vch\ 
of —0.91(7)% in comparison to the tree-level result of the total decay rate. For the exclu- 
sive B^ — )• vr"'' / (7) decays, we find a correction of 5.50(1.50)% to the total rate due to 
next-to-leading order effects. This translates into a reduction of |Vub| by —2.64(34)%. 

The leading uncertainties m. B ^ Dlv (7) and B ^ nliy (7) results originate from 
neglecting short-distance contributions and due to the matching procedure. The large 
matching uncertainty for B — t- 7r/i/(7) indicates, that the intermediate region between 
long- and short-distance regime is not modeled adequately with our simplified approach. 
The correction for iV^bl should therefore be used with care. This seems not to be the 
case for B ^ Dlv (7) , where swapping contributions between the long- and short-distance 
regime only provokes a modest scale dependence. Note that in addition to the estimated 
uncertainties both results might suffer from additional model dependencies, which cannot 
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be quantified within the model itself. The result for B — t- Dq I v (7) has the same feature, 
but due to the lack of further knowledge due to SD related real corrections, the pure IB 
result should be used with care. 

The extracted correction factors ioi B ^ Dlv (7) can be combined with the strong 
corrections to determine the corrections of the semileptonic decay rates at ©(a Gf+Os Gp). 
Corrections of mixing terms, which are at next-to-next-to-leading order ©(aa^Gp) are 
negligible. By providing the experimental community with these results and a next-to- 
leading order Monte Carlo generator, we hope to make a first step towards closing this 
important knowledge gap and improve the future extractions of | Vch \ and | Vah \ ■ 
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Table 11: The running of 1 + Ssd + <^id as a function of /io is shown ior B ^ Dlv (7). The averaged 
difference between the predictions of (Jtotai at ^/xo and 2/io to the prediction at /to = m£)+ or m^io is 
taken as the estimator due to matching uncertainty CTmatching- The uncertainty in the parentheses 
is numerical. 
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Table 12: The running of 1 + Ssd + <5id as a function of /io is shown ioi B ^ ttIi/ (7). The averaged 
difference between the predictions of ^totai at ^/io and 2/to to the prediction at /io = or m^o is 
taken as the estimator due to matching uncertainty (Tmatching- The uncertainty in the parentheses 
is numerical. 
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Table 13: The running of 1 + iJgd + ^id as a function of /xq is shown for for B — )• DqIu (7). The 
averaged difference between the predictions of c^totai at and 2^0 to the prediction at no = to^. + 
or m^i.o is taken as the estimator due to matching uncertainty Umatching- The uncertainty in the 
parentheses is numerical. 
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PI three- momentum: D'^ ei/ (j) pi three- momentmii: — ^ D'^ ^ f (7) 




0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 

/t [GeV] k [GeV] 



Figure 8: — > 13+ ^ (7): The predicted lepton and D ttiree-momentuni magnitudes, pi and p, 
for tree-level and next-to-leading order and the logarithmic photon energy k distribution are shown. 
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PI three- momentum: B'^ — > D^eiy{'i) pi three- momentum: B'^ — > D^jiul"/) 
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Figure 9: — > ^ (7): The predicted lepton and D three-momentum magnitudes, pi and p, 
for tree-level and next-to-leading order and the logarithmic photon energy k distribution are shown. 
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Figure 10: ^ -k^ Iv (7): The predicted lepton and tt three- momentum magnitudes, and p, 
for tree-level and next- to- leading order and the logarithmic photon energy k distribution are shown. 
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Figure 11: 5+ — > 7r°Zz^(7): The predicted lept- 
for tree-level and next-to-leading order and the loi 



PI three-momentum: B'^ — > tt^ pi^ (7) 




PI [GeV/c] 

p three-momentum: i?+ — ^ tt" (7) 




0.5 1 1.5 2 2.5 3 



p [GeV/c] 

k energy: £i+ — > tt'^ ji v 7 
T 1 =r 

■s - 
a : 

-H^lO-1 r- 



10-^ =- 




0.5 1 1.5 2 2.5 3 

k [GeV] 

in and tt three-momentum magnitudes, and p, 
;arithmic photon energy fc distribution are shown. 
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p; three-momentum: ^ D^'^ e u (j) pi three-momentum: ^ D'^^^ jiu (7) 
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Figure 12: — > + I v (7): The predicted lepton and Z?o three-momentum magnitudes, and 
p, for tree-level and next-to-leading order and the logarithmic photon energy fc distribution are 
shown. 
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Pi three-momentum: D^,^ t u (j) pi three-momentum: B+ — > D^'^ iiu{^) 
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Figure 13: — > Z)g° I v (7): The predicted lepton and Z?o three-momentum magnitudes, and 
p, for tree-level and next-to-leading order and the logarithmic photon energy fc distribution are 
shown. 
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A. Form factor models of exclusive semileptonic S-meson decays reviewed 



A.l Form factors for B ^ Dlv 



Considering the heavy quark sector of the Standard Model Lagrangian in the limit of 
infinitely massive quarks yields interesting simplifications: The arising effective field theory 



is known as heavy quark effective theory [35| and predicts the expansion of the hadronic 
current for exclusive B ^ Dlu decays as a function of the four- velocity transfer w = vb-vd, 
where vb and vd are the four velocities of the B- or L>-meson, respectively. The relevant 
vector current is 

= V^^^^(h+{w){vB + vnY + h.{w){vB-VDY), (A.l) 
with the heavy quark form factors 



h^{w) = g{l) X l-8pl,z + (51p|, - 10) z'^ - {252pl, - 84) z^' 
h^{w) = 0, 



(A.2) 
(A.3) 



where z = ^^=-j-^ . The expansion depends on pj^ and ^(1), which are the form factor 
slope and the normalization at = 1, respectively. The numerical value of the slope used 



in this analysis is listed in Table 14 



Parameter Value 
Pd 1-19 
Table 14: Used slope taken from p&. 



A.2 Form factors for B ^ D* Iv 

Similarly as for the pseudoscalar D-meson final state, HQET can predict the total and 
differential rates for B — )• D* I u decays as a function of the four- velocity transfer w = 
Vb • vd*- The vector and axial- vector currents are [35| 

{D*\V^\B) = 2ie^^^^e:{pB)SpD^)pg{t), (A.4) 
{D*\A^'\B) = e*^' f{t) + {pB+PD*Ye* ■pBa+{t) + {pB-PD*Ye* ■pBa.{t),{K.b) 



with 



= — inT—-\ ^1 ' (^-6) 

niB K[l + r) 

f{t) = — + 1) hAi (w) , (A.7) 

= a Mw) hAi{w) , (A.8) 

niB K[l + r) 

a-{t) = 0. (A.9) 
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and 

hAiiw)=T{l)x 1 - 8pl), z + {53pl, - 15) z"^ - {231pl), - 91) (A.IO) 

Riiw) = Ri{0) - 0.12{w - 1) + Om{w - 1)2 , (A.ll) 

R^{w) = i?2(0) + 0.11(?i; - 1) - OM{w - 1)2 , (A.12) 

where p'j^* and J-{1) are, respectively, the form factor slope and normalization at w = 1. 
The used numerical values of the slope and ratios are stated in Table 15. 

Parameter Value 

Pd 1-16 
i?i(0) 1.369 
i?2(0) 0.846 



Table 15: Used slope from 



A. 3 Form factors for B ^ ttIv 

Extrapolating results from lattice QCD calculations and light-cone sum rules, the form 
factors of exclusive B ^ ttIu can be parametrized as a function of the four-momentum 
transfer squared t. The relevant vector current is |37] 



(vriy^i?) = {{PB +p.r - "^^^j^iPB-p.r) hit) + { '^'\'^' {PB-p.)') Mt) , 

(A.13) 

with form factors parametrized as 

/+W = 7^^ + 7^^, (A.14) 

i T 1 T 

m, m, 
t+i J+2 

hit) = ^^^V ' (A-15) 

fo 

with the normalizations rf^-^,rf^^ , and rjg, and the pole masses m-y+i, rnj^^, and mj^^. The 
values used in this paper are stated in Table |l^. The proportional form factor f-{t) 
can be written as 



f-{t) = ^^^^j^{Mt)-U{t)) 



t 



~5 -L 2- 



m, nif m, 

to J+l J+2 



(A.16) 



A. 4 Form factors for B 

Heavy quark effective theory can be used to extract differential and total rates for other 
charmed resonances in an approximative manner, as done in 39 1 for B — )■ D^lv decays. 
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Parameter Value 





0.744 




-0.486 


2 

rUf 

j+i 


28.40 GeVVc^ 


2 

rrif 

i+2 


40.73 GeVVc^ 




0.258 




33.81GeVVc^ 



Table 16: Form factor slopes and pole masses taken from |37 



Since the Dq is a scalar particle, the hadronic vector current vanishes and only axial 
contributions are non-zero. Note that the entire formalism discussed in Sec. |3| can be 
adapted to derive the IB axial-vector contributions with SD related corrections. The 
relevant hadronic contribution is given by 



with the form factors g±. They are parametrized by 

wA* -A 



g+{w) 



2{w - l)Ci{w) - 3C{w) 



w + 1 

A* {2w + 1) - A{w + 2) 



w + 1 



CM-2(u;-l)CiH 



9-{w) = CM 



where 



CM 



w + 1 



t{w) . 



with t{w) is the leading Isgur-wise function pO| , |41[| . This can be simplified to 

CM = C(i)x [i + c'iw-i)]., 



(A.17) 



(A.18) 



(A.19) 



(A.20) 



where C' denotes the form factor slope. The value of Ci depends on the phase-space region, 
but its exact value only leads to small corrections. As a consequence, its contributions can 
either be neglected or approximatively taken into account. The latter is done by setting 

Ci = AC(^) • (A.21) 
The numerical values of the remaining slopes and pole masses are stated in Table 0. 



-46- 



Parameter Value 



Eb 0.1042 c GeV-i 

€c 0.3571 c GeV"^ 

A 0.4 GeV/c^ 

A* 0.75 GeV/c2 

C' -1.0 



Table 17: Form factor slopes and pole masses taken from |38 



B. Loop integrals 



B.l Passarino Veltman reduction of higher order tensor integrals 

For one-loop tensor integrals, a systematic algorithm has been worked out by Passarino 
and Veltman [27| to reduce any given higher order tensor integral as a sum of four- vectors 
multiplied by scalar integrals. For integrals with one, two, or three external legs, we use 
the notation 

A I 2\ ~A-D f k 1, fe^fcjy . , 



(2vr)^/2 (-A;2 + 



m 



D / 2 2 2\ ~4-D f '^^^ ^,kf^,kfj_ki, /-DON 

Bo,,MP ;-0,- ) = j (2,)D/2 (_,2 + ^2) ((_^ + ,)2+^2) ' (B-2) 



~4— D j d. k 1, kfj^, kfj^ki/ 



(2vr)^/2 (-fe2 + ml) {{-pi + kf + mf) ((-p2 + + ml) ' 

(B.3) 



respectively, with an apparent generalization towards more external legs and higher rank 
tensor integrals. Note that we redefined fi^~^ = —i (47r)^ and neglected A^. The 

latter vanishes due to its odd transformation property when integrated. In addition, we 
defined s = {pi — ^2)^- Due to Lorentz symmetry the integration result can only depend 
on tensor structures from external momenta p^ and the metric tensor g^^^ . This implies 

B'^- =p^^p-Bu+9^''Boo, 

= Pt p'i Cn + p'i C22 + P^^ C12 + 5^" C7oo . 



Contracting both sides of these equations with external momenta and the metric tensor 
yields a set of equations that allow the determination of the loop-form factors Sj, Bij,Ci, Cij 
with i,j = 0,1,2. 
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B.1.1 Derivatives of two-point functions 

The derivative of a two-point function is defined as 

Bi = ^Biip^;mlm^), (B.4) 

with i = 0,1. 
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